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Abstract.  Carbon  vapor  formed  at  graphite  evaporation  contains  considerable  amounts  of  molecules  Co  and  C3  along 
with  single  atoms  C.  A  numerical  kinetic  approach  to  gas  mixtures  is  proposed  which  is  based  on  a  relaxation  model  for 
the  Boltzmann  equation.  It  is  found  that  relations  between  the  pressure  and  temperature  ratios  and  the  Mach  number 
behind  the  Knudsen  layer  do  not  depend  on  vapor  composition  and  are  the  same  as  in  one-component  vapor.  The  molecu¬ 
lar  composition  of  vapor  is  determined  by  decreasing  the  thermal  velocity  with  molecular  mass  that  decreases  evapora¬ 
tion  rate  of  heavier  species  and  by  decreasing  the  recondensation  fraction  with  molecular  mass  that  acts  oppositely. 


INTRODUCTION 

Carbon  vapor  obtained  by  evaporation  of  graphite  is  known  to  contain  considerable  fractions  of  molecules  C2  and 
C  5  along  with  single  atoms  C  [1].  Recent  experiments  on  near-threshold  laser  ablation  [2]  indicate  that  the  molecules 
do  not  form  in  the  gas  phase  but  are  ejected  directly  from  the  condensed  phase.  The  molecules  may  act  as 
condensation  centers  that  is  very  important  in  gas-phase  synthesis  of  larger  clusters  [3],  fullerens,  and  single-wall 
carbon  nanotubes  [4].  Evaporation  rate  and  molecular  composition  of  vapor  both  depend  on  microscopic  laws  of 
molecule-condensed  phase  interaction  and  on  molecular  collisions  in  the  Knudsen  layer  in  vapor  adjacent  to  the 
surface.  The  present  work  is  devoted  to  a  multicomponent  Knudsen  layer. 

Strong  evaporation  problem  is  well  studied  for  one-component  vapors  [5-7]  .  Peculiarities  arising  at  evaporation 
of  two-component  substrates  have  been  analyzed  [7-8]  by  Direct  Simulation  Monte  Carlo  simulation.  An  approach 
to  binary  mixture  condensation  [9]  is  known  where  a  relaxation  model  for  gas  mixture  [10]  was  numerically  studied. 


MODEL 


Evaporation  from  a  plane  condensed-phase  surface  is  considered.  Relaxation  model  [10]  is  accepted  with  an  N- 
component  gas  mixture  described  by  velocity  distribution  functions  fa  satisfying  the  transport  equations: 


^  +  c=ff  =  2 >«p(^-/a)  for  a  =  1  ..  N, 


(1) 


M 


where  t  is  the  time,  z  the  coordinate  along  the  normal  to  the  surface,  cz  the  projection  of  the  molecular  velocity  to  the 
normal,  0)(,p  the  effective  collision  frequency  of  a  particle  of  species  a  with  particles  of  species  p.  and  Fap  the 
reference  distribution  functions.  The  collision  terms  satisfy  the  collisional  invariant  conditions: 

J(Ga-/a)CU;,C2)dc  =  0,  (2) 

J(Gp-/a)dc  =  0,  (3) 

J(Gp  -  fa)(cz,c2)&c  +  copamp  J(Fpa  -  /p)(c. ,c2)dc  =  0 ,  (4) 


where  c  is  the  vector  of  molecular  velocity,  c  its  absolute  value,  and  ma  the  molecular  mass  of  species  a.  Two 
additional  conditions  specific  for  Maxwell  molecules  are  required  [10]: 
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where  cr  is  the  projection  of  the  molecular  velocity  to  the  condensed  phase  surface,  k  the  Boltzmann  constant,  and 
na,  ua  and  Ta  are  the  number  density,  the  average  velocity,  and  the  temperature,  respectively,  of  species  a. 


Discrete  Velocity  Method 


Discrete  velocity  method  [11]  is  applied  for  numerical  solution  of  system  (1).  The  two-dimensional  velocity 
space  (c-,  c,)  is  divided  into  equal  rectangular  cells  with  center  coordinates  c\  and  c  ' ,  where  i  and  j  are  integers.  The 
discrete  velocity  distributions,  ,  are  defined  as  averages  of  fa  over  cell  (/,/).  Approximation  of  the  left-hand  side 
of  Eq.  (1)  by  a  second-order  essentially  non-oscillatory  scheme  is  described  in  [11]  while  evaluation  of  the  discrete 
reference  distributions  F‘L  is  original  in  the  present  work.  They  are  looking  for  in  the  following  form  corresponding 
to  the  continuous  model  [10]: 

F&  =  (Axp  +  -«)  +  CaliK<  -«)2  +  (c/)2]}exp[-ra[(ci  -  AJ2  +  (c/)2]} ,  (7) 


where  u  is  the  mass-averaged  velocity  of  the  mixture.  Coefficients  Ta  and  Aa  specify  equilibrium  discrete 
distributions  corresponding  to  molecular  masses,  ma ,  and  temperature,  T,  and  velocity,  u,  of  the  mixture.  They 
satisfy  the  following  equations  [11]: 


'XbJ 

R{=0,  R2S0  +  R0(S2  S0)  =  0, 

ma 

(8) 

with  the  sums 

RV=T  (<  ~  “)Y  eXP(-ra«  -  AJ2]  , 

(9) 

sy  =  ^(c  ’  )Ycxp|-ra(cj  )2  \Wj  , 

(10) 

j 


where  w;  ~  cJr  is  the  volume  of  velocity-space  cell  (i,  j)  [11].  Equations  (8)  are  numerically  solved  by  the  Newton 
method  with  the  initial  approximations  of  Ta  =  m,,/(2kT)  and  Aa  =  u,  which  correspond  to  the  continuous  limit. 

A  discrete  approximation  of  Eqs.  (2)-(6)  gives  the  following  system  of  linear  equations  for  each  triple,  Aap,  Bap, 
and  Cap,  at  any  combination  of  a  and  p: 

Axp^cA  ^aP^I^O  Cap  (F2S0  +  R0S2  )  =  1 1a  ,  (1 1) 

Axp^i^o  ^ap^2^o  CgpiRjSg  +  ^?15'2)  =  na(ua p  —  u)  ,  (12) 

Aap (^2^0  +  R<)S2 )  +  Bap (R}S0  +  RlS2)  +  Ca/i{R4S0  +  2R2S2  +  R0S 4)  =  3nakTap  / ma  +  na(uap  —u),  (13) 


where 


Vp 


maiia  +  mpKp 


+  m,. 


Taf,  -Ta  + 


2m„ma 


(ma  +  77/p  ) 


-[T’p-^  +  ^K-Up)2]. 


(14) 


Microscopic  Characteristics  of  Carbon  Vapor 

A  three-component  mixture  of  atoms  C  and  molecules  C2  and  C3  is  considered.  While  4-  and  5-atomic  molecules 
are  often  present  in  carbon  vapor  [1-2],  their  concentrations  are  negligible  and  so  they  are  not  taken  into  account. 
Collisions  between  molecules  may  result  in  changing  their  vibrational  and  rotational  states  and  even  in  their 
dissociation  or  merging.  Cross-sections  of  these  non-elastic  processes  are  supposed  to  be  considerably  less  than  the 
elastic  collision  cross-section.  Therefore,  non-elastic  collisions  can  be  neglected  within  the  Knudsen  layer  where  an 
evaporated  molecule  is  subjected  to  few  collisions  only.  Equation  (1)  actually  describes  translational  degrees  of 
freedom  only.  Internal  states  are  supposed  to  be  uniform  in  the  Knudsen  layer. 

To  estimate  relaxation  frequencies  0)„p,  a  correlation  [11]  between  the  rigid  spheres  model  and  the  relaxation 
model  is  used.  In  case  of  mixture,  the  relaxation  frequencies  are  proportional  to  the  corresponding  collision  cross- 
sections  aap  of  molecules  a  with  molecules  p: 
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(15) 


The  most  stable  configuration  of  molecule  C3  is  the  linear  one  [12],  therefore  it  is  considered  as  a  circular  cylinder 
with  hemispherical  ends  as  shown  in  Fig.  1  (a).  The  same  geometry  describes  molecules  C2  and  atoms  C  (reduced 
case  with  cylinder  length  equal  to  zero).  The  only  parameter  characterizing  size  of  the  molecules  is  atomic  diameter 
a.  This  value  is  taken  equal  to  2  Angstroms  that  roughly  corresponds  to  mean  interatomic  distance  in  graphite. 
Cross-sections  aKp  are  estimated  in  assumption  that  striking  molecules  can  be  oriented  along  the  three  coordinate 
axes  that  gives  9  possible  pair  combinations,  which  can  be  classified  as  shown  in  Fig.  1  (b).  Averaging  over  all  the 
orientations  gives  the  mean  collision  cross-section  of  a  molecule  Ca  with  a  molecule  Cp: 

aaB  =7t<72{l  +  —  [ap  +  5(a  +  p)-ll]}  for  a,p  =  1,  2,  3.  (16) 

97t 

Experimental  data  concerning  interaction  of  vapor  molecules  with  the  condensed  phase  surface  are  rather  limited 
[1,2].  Complete  sticking  of  impinging  molecules  is  supposed  in  the  present  model  that  corresponds  to  the 
condensation  coefficient  of  1.  Velocity  distributions  of  molecules  ejected  from  the  condensed  phase  are  specified  by 
Maxwellians  corresponding  to  surface  temperature  T,  and  partial  pressures  of  components  psa  in  the  saturated  vapor: 
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(17) 


RESULTS  AND  DISCUSSION 


Figures  2-3  show  detailed  calculation  results  for  an  example  of  steady-state  evaporation  at  the  pressure  ratio  pjps 
=  0.22  where  pa  is  the  ambient  pressure  far  from  the  Knudsen  layer.  Equal  partial  saturated  pressures  psi  =  ps2  =  ps 3  = 
pj 3  are  chosen  to  exhibit  peculiarities  caused  by  multicomponent  vapor.  A  strong  evaporation  flow  with  the  Mach 
number  M  —  0.964  at  infinity  is  formed  at  such  parameters. 

Figure  2  shows  profiles  of  characteristic  moments  of  velocity  distribution  functions  in  the  Knudsen  layer. 
Dimensionless  parameters  are  used  with  pressure  and  temperature  normalized  by  saturated  vapor  pressure  ps  and 
surface  temperature  Ts,  respectively.  Density  and  velocity  are  normalized  by  saturated  vapor  density  ns  =  ps/kTs  and 
thermal  velocity  of  monatomic  carbon  vapor  at  the  surface  temperature  u0  =  (2kTJm\)12 .  Net  heat  flow  q  in  the 
frame  moving  with  the  vapor  is  measured  in  the  units  of  «0 ps.  Dimensionless  coordinate  in  this  figure  is  the  ratio  of 
distance  z  from  the  surface  to  the  mean  free  path  of  monatomic  carbon  vapor  at  saturated  vapor  density  ls=  l/(na2ns). 


(a) 
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C2 


C3 


(b) 

cross-section  weight 

na2  1 

[jt  +  2(a  +  |3  -  2)]a2  2 

[7t  +  2(P  -  l)]a2  2 

[7t  +  2(a  -  l  )]a2  2 

[ji  +  (a  +  1)0  +  1)  -  4]a2  2 


FIGURE  1.  Geometric  parameters  of  carbon  molecules  (a)  accepted  for  estimation  collision  cross-sections  at  various  mutual 
orientations  of  striking  molecules  Ca  and  Cp  (b).  Symbol  a  designates  the  carbon  atom  diameter.  Collision  cross-sections  at 
particular  orientations  along  with  their  weights  are  indicated  on  the  right. 
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All  the  parameters  shown  in  Fig.  2  change  sharply  near  the  surface  but  tend  to  constant  values  at  infinity. 
Number  density  of  vapor  (see  broken  line  in  Fig.  2  (c))  and  all  its  components  (Fig.  2  (a))  decrease  in  the  Knudsen 
layer.  It  should  be  noted  that  number  densities  of  carbon  vapor  species  considerably  decrease  with  their  molecular 
masses  (see  Fig.  2  (a))  while  their  densities  in  the  saturated  vapor  are  equal  in  the  considered  case.  The  difference  in 
the  carbon  species  transport  in  the  Knudsen  layer  is  clearly  seen  in  Fig.  2  (b)  where  relative  concentrations 
xa  =  nj{n  i  +  n2  +  n3)  are  plotted.  The  concentration  of  atoms  is  higher  than  the  value  of  1/3  in  the  saturated  vapor 
and  increases  with  distance.  The  concentration  of  molecules  C2  is  about  1/3  and  that  of  molecules  C3  is  less  than  1/3 
and  decreases  with  distance.  This  can  be  qualitatively  explained  by  difference  in  thermal  velocities:  the  lightest 
monatomic  component  with  the  highest  thermal  velocity  evaporates  more  effectively  than  the  heaviest  molecule  C3 
with  the  lowest  thermal  velocity.  Similar  result  was  obtained  for  two-component  vapor  [7-8]. 

Parameters  of  the  multicomponent  vapor  in  the  Knudsen  layer  presented  in  Fig.  2  (c),  (d)  are  qualitatively  the 
same  as  for  one-component  vapor  [6,1 1].  Local  Mach  number  M  =  u/c  (upper  full  line  in  Fig.  2  (c))  is  defined  as  the 
ratio  of  vapor  velocity  u  to  local  sound  speed  c  =  | (5/3) pl(miii\  +  m2n2  +  m3n3)]12,  where  p  is  the  pressure  of  the 
mixture.  The  last  equation  for  the  sound  speed  is  obtained  in  assumption  of  no  exchange  between  translational  and 
internal  degrees  of  freedom  that  corresponds  to  the  accepted  kinetic  model  (1).  The  local  Mach  number  increases  in 
the  Knudsen  number  and  tends  to  the  value  of  0.964  at  infinity. 

Difference  between  axial.  pzz  (full  line  in  Fig.  2  (d)),  and  radial,  pn  (broken  line),  pressure  tensor  components  and 
deviation  of  heat  flow  q  (chain  line)  from  zero  characterize  degree  of  nonequilibrium  that  can  be  used  to  define 
thickness  of  the  Knudsen  layer  [11],  One  can  see  that  the  layer  of  “mechanical”  nonequilibrium  with  pzz=t  prr  is 
thinner  than  the  layer  of  “thermal”  nonequilibrium  with  q  ^  0  as  in  case  of  one  component  [11].  Knudsen  layer 
thickness  considerably  increases  when  Mach  number  at  infinity  approaches  unity  as  in  one  component  vapor  [11]. 

Figure  3  shows  velocity  distribution  functions  of  the  three  components  at  the  surface.  The  forward  part  of  the 
velocity  distributions  is  the  half-Maxwellian  defined  by  Eq.  (17)  and  the  backward  part  at  cz  <  0  characterizes  back 
condensation  of  evaporated  species  that  have  changed  their  directions  due  to  collisions  in  the  Knudsen  layer.  The 
backward  part  is  rather  sensitive  to  processes  in  the  Knudsen  layer  [5-6,11].  A  discontinuity  at  the  plane  cz  =  0 
between  the  forward  and  backward  parts  is  clearly  seen  in  Fig.  3  that  is  general  for  the  Knudsen  layer. 


FIGURE  2.  Distributions  of  dimensionless  parameters  in  the  Knudsen  layer  at  pressure  ratio  pjps  =  0.22  and  equal  saturated 
pressures  of  the  components  at  the  surface,  ps \  =  ps2  =  ps 3  =  pj 3:  (a),  dimensionless  densities  njns  of  the  components;  (b),  relative 
concentrations  xa  of  the  components;  (c),  local  Mach  number  M  and  dimensionless  density,  nlns,  velocity,  u/u0,  and  temperature, 
777),  of  the  mixture;  (d),  dimensionless  components  of  the  pressure  tensor,  pnJps  and  pzz/ps ,  and  heat  flow  q/(u0ps). 
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FIGURE  3.  Velocity  distributions  (arbitrary  units)  of  carbon  species,  C,  Co,  and  C3,  on  the  condensed  phase  surface  at  pressure 
ratio  pjps  =  0.22,  Mach  number  M  =  0.964,  and  equal  saturated  pressures  of  the  components.  Axial,  c.,  and  radial,  cr,  molecular 
velocities  are  normalized  with  thermal  velocity  u0  =  (2kT Jm\)m .  Recondensation  fraction  is  shown  on  each  diagram. 


Importance  of  the  backward  part  of  velocity  distribution  relative  to  the  forward  one  decreases  with  molecular 
mass  (see  Fig.  3).  This  is  in  line  with  molecular  collision  mechanics:  heavier  molecules  are  rarely  scattered  to 
backward  directions.  The  recondensation  fractions  of  carbon  species  shown  in  Fig.  3  are  calculated  as  ratios  of 
backward  molecular  flux  at  cz  <  0  to  forward  flux  at  cz  >  0.  Notice  that  the  recondensation  fraction  in  monatomic 
vapor  at  similar  conditions,  M  =  1,  is  0.185  [5],  Thus,  fraction  of  atomic  carbon  backscattered  in  the  Knudsen  layer 
is  higher  than  in  one-component  case  and  fractions  of  backscattered  molecules  C2  and  C3  are  less.  The  same 
tendency  of  increasing  the  recondensation  fraction  with  molecular  mass  was  observed  for  binary  mixture  [8], 

Figure  4  shows  calculations  at  different  pressure  ratios,  pjps,  between  0.22  and  1  that  almost  cover  the  entire 
range  of  Mach  numbers  from  0  to  1.  The  vapor  characteristics  in  this  figure  are  the  same  as  in  Fig.  2  but  they  refer  to 
the  region  behind  the  Knudsen  layer  where  there  is  no  dependence  on  coordinate.  Vapor  temperature,  Ta,  and  total 
number  density,  na,  and  pressure,  pa,  of  all  the  components  are  compared  with  the  corresponding  values  for  one- 
component  vapor  in  Fig.  4  (a).  It  is  noticeable  that  present  multicomponent  model  (crosses)  gives  the  same  results  as 
previous  calculations  (open  and  dark  symbols)  for  the  relaxation  model  in  one-component  case  [6,11].  The 
difference  is  within  numerical  error.  Similar  conclusion  was  obtained  for  binary  mixture  [8], 


FIGURE  4.  Carbon  vapor  parameters  far  from  the  surface  as  functions  of  the  Mach  number  at  equal  partial  pressures  p,t  =  ps2  = 
Psh  -  Ps! 3  in  the  saturated  vapor,  (a),  comparison  of  temperature,  TJTS,  density,  «„/«.„  and  pressure,  p„/ps,  ratios  for 
multicomponent  carbon  vapor  (crosses)  with  those  for  one-component  vapor:  curves,  analytical  theory  by  Ytrehus  [5];  open  and 
dark  symbols,  numerical  calculations  for  the  relaxation  model  by  Sone  et  al.  [6]  and  by  Gusarov  et  al.  [11],  respectively,  (b), 
relative  concentrations  of  carbon  species,  xaa:  symbols,  numerical  results;  full  line  curves,  their  spline  interpolations;  broken 
lines,  collisionless  evaporation  (the  Hertz-Knudsen  equation). 
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Essential  results  concern  molecular  composition  of  carbon  vapor  presented  in  Fig.  4  (b).  At  low  Mach  numbers, 
M  «  1,  vapor  composition  is  about  the  composition  of  the  saturated  vapor  with  xal  =  xa2  =  xa2  =1/3.  Concentration 
of  carbon  atoms,  xal,  increases  with  the  Mach  number  because  of  their  higher  thermal  velocity  relative  to  the 
molecules.  However,  xal  does  not  attain  the  collisionless  limit  given  by  the  Hertz-Knudsen  evaporation  low  (upper 
broken  line  in  Fig.  4  (b)).  This  appears  to  be  due  to  the  highest  recondensation  fraction  of  the  atoms.  On  the 
contrary,  concentrations  of  carbon  molecules,  C2  and  C3,  decrease  with  the  Mach  number  tending  to  their 
collisionless  limits  (middle  and  lower  broken  lines)  but  do  not  attain  them. 

Thus,  the  molecular  composition  of  carbon  vapor  is  determined  by  two  opposite  tendencies:  the  first  is 
decreasing  the  thermal  velocity  with  molecular  mass  that  leads  to  decreasing  the  concentration  and  the  second  is 
decreasing  the  recondensation  fraction  with  molecular  mass  that  leads  to  increasing  the  concentration.  Actually,  the 
first  tendency  is  more  important  for  the  considered  conditions,  and  so  the  concentration  of  carbon  species  decreases 
with  their  molecular  mass. 

The  considered  example  of  graphite  evaporation  at  equal  partial  saturated  pressures  of  C,  C2,  and  C3  suggests  a 
general  conclusion  that  evaporation  rate  and  vapor  temperature  do  not  depend  on  molecular  composition  of  vapor 
and  can  be  obtained  from  the  well-known  theory  of  one-component  evaporation.  On  the  other  hand,  the  molecular 
composition  itself  may  considerably  differ  from  the  composition  of  the  saturated  vapor  and  can  be  obtained  using 
the  present  multicomponent  kinetic  model. 


CONSLUSIONS 

A  numerical  kinetic  approach  to  gas  mixtures  is  proposed  which  is  based  on  a  relaxation  model  for  the 
Boltzmann  equation.  Strong  evaporation  of  graphite  with  formation  of  a  three-component  vapor  is  analyzed  at  equal 
saturation  pressures  of  atoms  C  and  molecules  C2  and  C3. 

The  structure  of  a  multicomponent  Knudsen  layer  at  strong  evaporation  is  qualitatively  the  same  as  in  the  one- 
component  case.  However,  molecular  composition  of  the  vapor  changes  through  the  Knudsen  layer  when  the 
distance  to  the  surface  increases,  and  tends  to  a  composition  different  from  that  of  the  saturated  vapor. 

Relations  between  the  pressure  and  temperature  ratios  and  the  Mach  number  behind  the  Knudsen  layer  do  not 
depend  on  vapor  composition  and  are  the  same  as  in  one-component  vapor.  The  molecular  composition  of  vapor  is 
determined  by  decreasing  the  thermal  velocity  with  molecular  mass  that  decreases  evaporation  rate  of  heavier 
species  and  by  decreasing  the  recondensation  fraction  with  molecular  mass  that  acts  oppositely. 
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